A dynamical model of the chemical bond: 
kinetic energy resonances between atomic orbitals 
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A new approach to chemical bonding is introduced in 
order to provide an improved understanding of the connec- 
tion between basic quantum mechanics and the covalent pair 
bond. It's focus is on the fact that the energy of the bond 
is largely given by the kinetic energy of the electrons, while 
the Coulomb forces are only comparable to the kinetic en- 
ergy terms close to the atomic nuclei, where they define the 
shape and the size of the atomic orbitals. If atomic orbitals 
are used as a starting point, the kinetic energy operator is 
sufficient to determine the energy of the chemical bond. The 
simple mathematical structure of this operator allows a cal- 
culation of bond energy as a function of the distance between 
the atoms. For Gaussian wavefunctions, it is possible to calcu- 
late this bonding potential analytically, determining the bond 
length, the bond energy and the elasticity of the bond from 
only a single parameter, the width of the atomic wavefunc- 
tion. It is shown that the results correspond surprisingly well 
with experimental values for diatomic molecules. 



I. INTRODUCTION 

While the theory of chemical bonding has advanced 
rapidly as more and more details of the quantum mechan- 
ical problem have been taken into account, surprisingly 
little effort has gone into the developement of intuitive 
models to interpret the results obtained by more sophis- 
ticated means. As a result, the theory of electronic struc- 
ture is still introduced quite differently to chemists (for 
example jlj) and to physisists (for example Q). While 
the former begin with energetic considerations in small 
molecules, the latter tend to focus on dynamics and solid 
state systems. Yet, the nature of the bonds is quite sim- 
ilar in solids and in molecules, and a lot can be learned 
by examining common features. 

The Sommerfeld theory of electrons in metals started 
out by completely ignoring the Coulomb interactions of 
electrons and nuclei, treating the electrons as free parti- 
cles moving in an effective vacuum. At first sight, this 
free electron nature of the metallic bond seems to be qual- 
itatively different from covalent bonds. However, in 1979 
Froyen and Harrison |3| showed, that the parameters for 
the LCAO or tight binding approach to the band struc- 
ture of tetrahedrally bound solids such as Diamond, Sil- 
icon and Germanium can be derived from a comparison 
with free electrons. Consequently, the matrix elements of 



the bonding overlap depend only on the kinetic energy. 
The potential energy enters through the separation of the 
atomic energy levels and through the bond length, both 
of which are features of the atomic wavefunction. 

It is therefore possible to explain bonding between 
atoms using only the contribution of the kinetic energy 
to the bond formation, while the Coulomb energy is re- 
sponsible for tightly binding the electrons to the atoms. 
This qualitative difference between the energy of atomic 
orbitals and the energy of molecular bonds is crucial for 
understanding the physics of the chemical bond. 



II. MECHANICAL PROPERTIES OF SINGLE 
ATOMS: A SIMPLE APPROACH TO 
UNDERSTANDING COMPRESSIBILITY 

Electrons are fermions. They obey the Pauli princi- 
ple. Therefore, two electrons of the same spin can only 
occupy the same region of space if they have different 
momentum. As a result, an increase in electron density 
requires an increase in kinetic energy and therefore elec- 
trons resist compression even if we neglect the Coulomb 
repulsion. In fact, for high densities, the Pauli principle 
turns out to be a stronger force of repulsion than the 
Coulomb forces. 

The correct quantitative expression can be derived 
from the Sommerfeld quantization rule that the density 
of quantum states in phase space is h~ 3 . The electron 
spin doubles this density of states, as there are two quan- 
tized spin directions for each phase space volume of h 3 . 
The minimal energy for a given density of electrons may 
now be determined by filling the lowest energy regions 
of phase space and calculating the average kinetic energy 
Epauii of each electron as a function of the electron den- 
sity p. The result, usually derived as part of the Sommer- 
feld theory of electrons in metals (for example in chapter 
2 of U), is: 
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where m is the electron mass and p is the electron 
density. 

This term, which increases with decreasing distance as 
r~ 2 , prevents the Coulomb attraction, which increases 
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only with r _1 , from pulling the electrons into the atomic 
nuclei. 

The details of an atomic system may be quite compli- 
cated. The exact density distributions of the electrons 
are subject to both quantum mechanical interference ef- 
fects and complicated screening effects which result from 
the electron-electron interaction. However, the basic me- 
chanical properties of size and compressibility may be 
explained without examining these details, just as it is 
done in macroscopic systems. For this purpose, we need 
to define a kind of mean field potential for the Coulomb 
interactions, including electron- electron repulsion. 

Effectively, there is one positive charge in the nu- 
cleus for each electron. Therefore, we will estimate the 
Coulomb energy of a given electron density by assuming 
that the charge of each electron is evenly distributed in a 
sphere with volume p^ 1 around a single positive charge. 
The Coulomb energy per electron is then given by 

Ecoulomb — ~Z ("^"P) 1 ^ 3 ' (2) 

87reo o 

Note that by ignoring the electrons which are much 
closer to the nucleus than most others and therefore see 
the "naked" charge of the nucleus we have essentially 
chosen to describe only the valence electrons. In fact, 
the "nucleus" is effectively an ion with a closed shell. 
The properties of size and elasticity are not affected by 
this ion, however, and we can therefore ignore this hidden 
part of the electronic system in this context. The choice 
of only a single positive charge for each electron should 
be best for alkali metals, so quantitative results should 
be tested by comparison with these elements. 

The complete thermodynamic relation for the energy 
E e i per electron as a function of the volume V e i = p^ 1 
per electron is 

P (V \ ~ 3k2 ( 3n2 \2/3 3e2 i 471 \l/3 / o\ 

10m Vei 8ire 3V e i 

The volume occupied by each electron in equilibrium 
can be determined by finding the minimum of the energy 
E e i. This minimum is found by setting the derivative 
dE/dV to 0. Note, that the volume derivative of the en- 
ergy is a pressure. Effectively, we determine the point of 
equilibrium between the Pauli pressure and the Coulomb 
pressure. In the atom, the electrons are compressed by 
the Coulomb attraction until the repulsion caused by the 
Pauli principle compensates this force. 

The length scale at which this happens is, of course, 
the atomic length scale of the Bohr radius a®. 

ao = l^!2j! = 0.529 x 10- 10 m (4) 
me 

The volume per electron calculated from the relation 
given above is V e i = 2(37r 3 ao/5) 3 w 47raQ, roughly equal 
to the volume of a sphere with radius 1.5ao- However, 
the approximation of the Coulomb pressure used is quite 
rough, so that only the order of magnitude is meaningful. 



To get an idea of the pressures involved, we can now 
calculate the Pauli pressure and the Coulomb pressure 
separately to find that for V e i ~ 47rap, this preassure is 
around 10 1 2Pa or 10 million times atmospheric pressure. 

To get an idea of how "solid" the atom actually is, we 
can estimate the compressibility K of the atom, defined 
as 

Not surprisingly, the result is roughly equal to the 
inverse Pauli pressure. For V e ; « 47raQ, we find K s» 
5 x 10~ 12 Pa~ 1 . This is exactly the order of magnitude 
found in the compressibilities of the alkali metals. For 
metals in general, typical values range from caesium at 
7 x lO^Pa- 1 to copper at 0.0745 x KT^PeT 1 . Note 
that the compressibility of metals does not depend on 
the properties of the metallic bond: we are not consider- 
ing the effects which make atoms stick together yet. To 
check, whether atoms which form only weak bonds show 
the same mechanical properties as the atoms of metals, 
a comparison with the noble gases is useful. The solids 
formed by condensing the noble gases at very low tem- 
peratures have compressibilities from 9 x 10 - 12 Pa -1 for 
neon to 2.8 x 10~ 12 Pa -1 for xenon. This is indeed very 
close to the alkali metals. Another type of crystal in 
which the hard sphere model of atoms seems to be ap- 
propriate is the ionic crystal. Although the charge dis- 
tribution in those crystals is very inhomogeneous, their 
compressibility is only one order of magnitude less than 
that of the alkali metals and the noble gases, with typical 
values close to the compressibility of 0.42 x 10 _12 Pa _1 
for sodium chloride. All values are taken from the corre- 
sponding chapters of 0. 

The result that the compressibility of single atoms is 
close to the compressibility of solids is well in tune with 
the quite intuitive picture of solids consisting of densly 
packed atoms. As a result of the immense magnitude of 
the Pauli pressure and the Coulomb pressure, the hard 
sphere model is quite an adequate picture for the atom 
despite the misleading notion of emptiness associated 
with the assumed "sizes" of electrons and nuclei. We 
can therefore conclude that the source of the low com- 
pressibility of solids is indeed the Pauli pressure and that 
the atomic properties defined by the equilibrium between 
Coulomb attraction and Pauli pressure are not changed 
very much when chemical bonds are formed. 

However, the Pauli pressure and the Coulomb pressure 
may be compensated atom by atom, and nothing in this 
picture makes the seperate hard sphere atoms stick to- 
gether. This is approximately true for the noble gases, 
but certainly not for most atoms. The periodic table of 
the elements shows us that the reacivity of the melements 
depends on the exact number of electrons. Therefore, we 
need to take a closer look at the particle properties of the 
electronic system in order to understand the forces that 
hold atoms together. 
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III. UNPAIRED ELECTRONS 

Since electrons are particles with spin 1/2 and can 
therefore occupy two spin states for every real space state, 
the Pauli pressure seems to act seperately on two types 
of electrons. If the electron gas were continuous and did 
not consist of quantized particles, the spin up density 
would thus always equal the spin down density, and all 
matter would be paramagnetic. Indeed, this type of para- 
magnetism exists in metals and is consequently known as 
Pauli paramagnetism (^]. 

However, in atoms there are a discrete number of elec- 
trons in a discrete number of states. Therefore, it may 
happen that n degenerate states are filled with less than 
2n electrons. In this case, the Pauli preassure is in- 
sensitive to the distribution of the electrons. Indeed, if 
all unpaired electrons align their spins, the Pauli princi- 
ple keeps them spatially separate, reducing the electron- 
electron Coulomb repulsion, while electrons with oppo- 
site spins would overlap more strongly, even if in separate 
real space states. This energy is the spin exchange en- 
ergy calculated in Hartree-Fock theory. If we have a high 
degeneracy of states partially filled, we therefore have a 
ground state with all unpaired spins parallel. This is the 
reason for both Hund's first rule and for ferromagnetism. 

Energetically, the unpaired electrons are always some- 
what unstable, as are all systems with a degenerate 
ground state, since any small perturbation will remove 
the degeneracy and break the symmetry of the system 
by forming new geometric patterns. In surface science, 
this tendency for forming bonds is often referred to as 
Jahn- Teller instabitity 0. Although the change in ge- 
ometry is not as impressive in the liquid or gas phase, 
the reactivity of atoms with partially filled shells is a 
consequence of an instability of the same type. 

The effect which causes the perturbation when an- 
other similar atom approaches is not, however, the elec- 
tric field, although van-der-Waals forces will be present. 
Instead, homopolar bond formation is largely the effect of 
kinetic energy, as the nearly free electron bands in solid 
state physics, which are found even in non-metallic sys- 
tems, suggest. The very small effects of charge density 
redistribution in a molecule are quite insufficient to ex- 
plain the energy of chemical bonds, especially since the 
energetically favorable positions are always the positions 
close to the atomic nucleus and not the position between 
the nuclei. In textbooks, this is often explained in a mis- 
leading way, because it seems so much harder to visualize 
the effects of kinetic energy in a picture. If the proper- 
ties of atoms that the kinetic energy and the potential 
energy are related by the equilibrium of Pauli pressure 
and Coulomb pressure are applied to the equations, it is 
even possible to express the kinetic energy contribution 
to bonding in terms of the potential energy of the atomic 
electrons. To avoid this type of confusion, We will there- 
fore neglect van-der-Waals forces and polar contribution 
altogether and determine the effect of kinetic energy only. 



This can also be justified by examining the Hamilto- 
nian of the two atom system: 

H = f + Vi + V 2 + V el - el (6) 

T is the kinetic energy of all the electrons, V1/2 is the 

potential of the respective atomic nucleus, and V e /_ e / 
is the electron-electron interaction. Considering that 
due to a negligible overlap of the electron distribution 
the Coulomb energies, including electron-electron inter- 
actions, act almost exclusively on wavefunctions localized 
at the respective atoms, this Hamiltonian may be seper- 
ated into atomic Hamiltonians: 

H 1 =f + V 1 + V eL _^ 1 (7) 

H 2 = f + % + V el - ell2 - (8) 

The eigenstates of Hi and Hi are the electronic states of 
the individual atoms. However, only the potential energy 
terms are localized at the site of the atom. The kinetic 
energy term T is the sum of all electronic kinetic energies 
regardless of the atom at which the electron is located. 
The total Hamiltonian is therefore given by the sum of 
the atomic Hamiltonians minus once the kinetic energy 
term. 

H = Hi+H 2 -f (9) 

For electronic states located entirely at atom 1, the con- 
tributions of H 2 — T are negligibly small and so are con- 
tributions of Hi — T for states located entirely at atom 
2. For delocalized electronic states however, equation (||) 
is a useful expression for the evaluation of energy contri- 
butions due to the delocalization as we shall see below. 

Of course, we have thus neglected the electron-electron 
interactions between electrons located at different atoms, 
as well as interactions between an atomic nucleus and 
electrons located at the other nucleus. For well sepa- 
rated atoms, the net effect of these interactions are the 
van-der-Waals forces and ionic interactions. We have ne- 
glected these effects, because van-der-Waals forces are 
known to be much weaker than the forces involved in 
chemical bonding and ionic effects are expected to be 
weak in the case of symmetric homopolar bonds. To 
treat the case of polar bonds, the ionic interaction term 
must be considered in more detail. 

Many particle effects within the atoms are still in- 
cluded in the Hamiltonian. However, they do not con- 
tribute to the bond energy. The kinetic energy term re- 
sponsible for the bond formation is generated because the 
kinetic energy is included in both atomic Hamiltonians, 
and must therefore be subtracted once from the total. 
This indicates the fact, that the kinetic energy is not lo- 
calized at one of the atoms, but causes a delocalization 
connecting the two atoms. 

The kinetic energy is a single particle property. There- 
fore, we can continue by calculating the energy contri- 
butions for each electron separately, adding the energy 
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terms of all electrons involved in the bond formation to 
obtain the total energy. The single particle wavefunctions 
may be obtained by Hartree-Fock theory, together with 
the effective single particle Hamiltonian for each atom. 

The matrix element connecting two electronic orbitals 
at different atoms, |1 > and |2 >, is given by 



< \\H e M\2 >= 
< l\Hf' f \2 > 



< l\HZ ff \2 > - < l\f\2 >, 



(10) 



where is the effective single electron Hamiltonian of 
atom i. Since |1 > and |2 > must be eigenstates of the re- 
spective atomic Hamiltonians the first two contributions 
are just the eigenenergies multiplied with the direct over- 
lap < 1|2 >. The latter is a correction for the fact that, 
since |1 > and |2 > are not orthogonal, a part of the 
matrix element represents the expectation value of the 
Energy of a state localized at one of the two atoms. If 
the non-orthogonality is removed, this contribution goes 
to 0, and the remaining matrix element is purely kinetic: 



< l\H eff \2 > 



< 1ITI2 > 



(11) 



Note, that the sign already suggests a lower energy for 
the symmetric state! 

It is now possible to develope a complete theory of 
bonding by defining the atomic wavefunctions and their 
energies and calculating the tight binding matrix ele- 
ments using only the operator of the kinetic energy. 



IV. BONDING AND ANTI-BONDING STATES 

If we consider only one electronic orbital per atom, the 
bonding and anti-bonding states of a homopolar bond 
must be the sum and the difference of the atomic wave- 
functions. The energies of the bonding and anti-bonding 
states can then be determined by calculating the expec- 
tation values of the kinetic energy for these two linear 
combinations. Note, that the sum and the difference of 
two normalized wavefunctions are automatically orthog- 
onal to each other. If the atomic wavefunctions are given 
by 1P1/2 (r), 
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= (V>i(r + d/2)±V; 2 (r-d/2)) (12a) 



^±(d) 



N ± = 2 ± /(Vi(r + d/2)V 2 *(r - d/2) 

+ Vi*(r + d/2)^ 2 (r-d/2))d 3 r 

1 



(12b) 



< Ekin >± (d) — — — (< Ekin >1 + < Ekin >2 



T 2^ I ^d(Mr + d/2))grad(r 2 (r - d/2)) 



+ grad(iPl(r + d/2))grad(ip 2 {r - d/2)))d 3 r (13) 

If each atom contributes one unpaired electron, these 
two electrons can redistribute into the two bonding and 
two anti-bonding levels as the atoms approach each other. 
If they pair up in the bonding level, a pair bond is formed. 

Although these real space integrals already describe 
the effect of wavefunction overlap on the kinetic energy 
completely, it is quite revealing to take a look at the same 
equations in k space by Fourier transforming the whole 
integral, which represents a convolution of two gradients. 
In k space, this changes into the integral of the product 
of the wavefunctions multiplied by k 2 and a cosine which 
represents the effect of the real space separation: 

< E km >± (d) = 



< Ekin >2 



T 2^ / (k 2 0Mk)V> 2 *(k)e- 



-ikd 



-^(k)V> 2 (6/fc)e+^ fcd ))d 3 k 



(14) 



If ipi = '02, the kinetic energy contributions of the 
overlapping wavefunctions can be written in an even more 
compact form: 

< E kin >± (d) = 



jL^L |(2k 2 0(k)^(k)(l± cos (kd)))d 3 k 



1% 2 /(k 2 ^(k)0*(k)(l ±cos(kd)))d 3 k 
2^ /(^(k)0*(k)(l±cos(kd)))d 3 k 



(15) 



In this case, the kinetic energy distribution is there- 
fore modified by a factor of 1 ± cos (kd) at each point 
in k space. If the width of the impulse distribution of 
the wavefunctions, that is, the width in k space of 
is roughly equal to ir/d, then the contributions with the 
highest k 2 values at the edge of the distribution are sup- 
pressed by a multiplication with values close to for the 
bonding state. The kinetic energy of the bond is therefore 
at a minimum. If the k space width of the atomic wave- 
functions is much larger than ir/d, the rapid oscillations 
of the cosine make all overlap contributions cancel (as 
they should, since the atoms are far apart in real space 
now). If the k space width is much smaller than ir/d, 
all parts of the k-space distribution contribute equally 
and the bonding state is again equivalent to the atomic 
wavefunction B. 

The k space distribution of the atomic wavefunction 
therefore defines a bond length of roughly ir divided by 
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its width in k-space, at which the average wavelength of 
the electron is in resonance with the bond length. Since 
the k-space width is connected with the real space width 
by the uncertainty relation, this bond length is roughly 
equal to the sum of the atomic radii, so the covalent bond 
is actually strongest when the atoms just about touch 
each other. 

To understand the significance of this result, it should 
be remembered, that the potential energy effects only 
enter the picture indirectly, by forming the atomic wave- 
functions. The formation of the chemical bond can then 
be explained entirely by the change of kinetic energy 
when unpaired electrons begin to tunnel resonantly be- 
tween the atoms. While the total electron density in real 
space changes only little, the kinetic energy distribution 
looses a major part of its high energy contributions. 



V. A QUANTITATIVE EXAMPLE: GAUSSIAN 
WAVEFUNCTIONS 

As an analytically solvable example, Gaussian wave- 
functions offer a good insight into the formation of bonds 
by kinetic energy resonance. The most simple bond is a 
symmetric combination of s-type states: 



V>(r) 



1 



(27T(7 2 ) 3 /4 



-r 2 /4<7 2 



(16a) 



Since the Fouriertransform of the Gaussian is again a 
Gaussian, it is easily possible to determine all features 
of the resulting bond between two unpaired electrons in 
such atomic states. 



^(k) = ( ^)3/4 e -^ 
7T 



< Ekin > — 



3h 2 

Ama 2 



(16b) 



(16c) 



Note that the factor of three is a consequence of 
the three spatial dimensions contributing to the energy. 
Since we can separate the dimensions for Gaussian wave- 
functions, we could as well examine only the one dimen- 
sional problem. However, we shall include the constant 
energy contributions for the sake of completeness. 



,-l/2{d/2a) 2 



l(dk/2 )(^!)3/4 f 
7T 



(18a) 



< Ekin > — 



37T 



3 - (3 - (d/2a) 2 )exp(-l/2(d/2a) 2 ) 



Ama 2 1 - exp(-l/2(d/2a) 2 ) 



{d/2af 



Ama 2 



Ama 2 exp{+l/2(d/2a) 2 ) - 1 



(18b) 



Figure 1 shows the energy and the k-space distribution 
of the bonding state as a function of d and figure 2 shows 
the same for the anti-bonding state. 
We can now determine the bond length, as well as the 
harmonic part of the potential around it. This allows an 
estimate of typical molecular vibrations and a compari- 
son with experimental values. 

The minimum of the bonding potential is at d — 3.2a. 
This corresponds very well with the estimate given in the 
previous section, that the bond length should be roughly 
equal to 7r divided by the width of the distribution in k- 
space. The standard deviation of the k-space distribution 
is l/2cr, so l/<7 is a good measure of its width. 
Around the minimum at d = 3.2a, the bonding potential 
may be written as 



< Ekin > — 



r((3 - 0.557) + 0.55(d/2<j - 1.6) 2 ). 



(19) 



The total energy is then given by Z times this potential, 
where Z is the number of electrons in the bond. From this 
potential, the following relations for bondlength do , bond 
energy Eb and the bond elasticity k can be obtained: 



dn = 3.2a 



E b = 51( 
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4.2( 



10 
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(20a) 



(20b) 



(20c) 



^+(k) 



-l/2(d/2 C r) 2 



C0S (dk/2)( — f'\ 
7T 



-<x 2 k 2 



< Ekin > z 



Ama 2 



3h 2 
Ama 2 



(17a) 



3 + (3 - {d/2a) 2 )exp{-l/2{d/2af 



l + exp(-l/2(d/2a) 2 ) 



{d/2af 



Ama 2 exp{+l/2(d/2a) 2 ) + 1 



(17b) 



Table 1 lists a few examples of diatomic molecules with 
their bond length, bonding energy and the elastic con- 
stant. Data for the elastic constant was taken from || 
and all other data from |J. The table also shows the 
a values that would correspond to the respective bond 
properties if the wavefunctions were simple Gaussians. 
Although there are strong deviations even between dif- 
ferent properties of the same molecule, the order of mag- 
nitude is reproduced correctly despite the simplicity of 
the model. 

To simulate a directional bond, we can calculate the 
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bonding potential for two wavefunctions of p-type sym- 
metry: 



VI. THE IONIC CONTRIBUTION: EXTENDING 
THE MODEL TO POLAR BONDS 



V>(r) 



^(k) 



0-(27TCT 2 ) 3 / 4 



2 I a 2 

-r /4cr 



2ak x ( — ) 3/4 e 

7T 



<E, 



kin 



>= 



5h 2 

Ama 2 



(21a) 



(21b) 



(21c) 



Note, that the separation between the two atoms, d, 
need not be along the x axis. It is possible to calculate 
the potential for the additive and the subtractive linear 
combinations as before. The potential for the bonding 
state in three dimensions is given by 



bh 2 

Ama 2 



h 2 (d x /2a) 2 ((d/2a) 2 -3 
Ama 2 4 - (d/2f7) 2 - exp(+l/2(d/2a) 2 ) 



(d y /2a) 2 + (d z /2a) 2 
exp(+l/2{d/2a) 2 ) - 1 



(22) 



This equation describes a potential with two minima 
along the x axis, at d x = ±5. 04c, as is shown in figure 3. 
It represents the bonding potential of a ppa bond, since 
there is rotational symmetry around the x-axis. In the 
vicinity of these minima, the potential may be approxi- 
mated by a harmonic potential of 

fi 2 

< Ekin > z 



Ama 2 



(5 - 0.73 + 1.08(-^ - 2.52) 
2a 



°.42((^) 2 + (^) 2 )). 

(23) 



It is now possible to compare the longitudinal and the 
transversal elasticity of the bond. The ratio of the cor- 
responding coefficients is 1.08/0.42. The bond is there- 
fore softer against shear forces than against compression 
by a factor of about 2.5. This example shows how the 
geometrical structure of molecules and their vibrations 
can be described using only the kinetic energy and the 
symmetry of atomic wavefunctions. Of course, a more 
comprehensive description of molecules or crystals would 
require the full calculation of matrix elements between all 
the atomic wavefunctions involved, using more realistic 
wavefunctions than the Gaussians presented here. 



The kinetic energy only dominates the homopolar 
bond. In the case of heteropolar bonds, there will be an 
energy difference between the electronic states at atom 
A and at atom B, which pulls the electron towards the 
more electronegative one. In a very simpleminded tight 
binding approach, this effect may be included as a site 
dependend energy. With Ea — Eb = 2D, the two by two 
matrix of the polar bond is then given by 



T+ D 
D T_ 



(24) 



T± is the kinetic energy of the bonding and anti-bonding 
states. For large separations d, this is an unrealistic ap- 
proach, however, since it neglects the energy needed to 
ionize the atoms and therefore predicts complete ioniza- 
tion as d approaches infinity, not including the Coulomb 
attraction. 

The other extreme is to assume fully polarized atoms. In 
this case, a term representing the Pauli repulsion is neces- 
sary to keep the ions apart. This may be done by adding 
the bonding and anti-bonding energies for the filled outer 
shells of the ions. 

To combine covalent and ionic effects into a model of the 
polar bond, it is necessary to consider a total of three 
possible electronic configurations: state 1, with both elec- 
trons at atom A, state 2 with one electron at each atom, 
and state 3 with both electrons at atom B. Note that 
this corresponds to a simple extension of Heitler-London 
theory ||. The Hamiltonian matrix is 

T/2{T+ - T_) 

T++T_ y/Tj2{T + -T-) 
1/2(T+-T_) T+ + T- + V B 

(25) 

For Va = +D and Vb = —D, this is equivalent to the two 
state model above. However, by including the ionization 
energy in the Va/b, we can now correct the result for 
large d. For very large d, Va/b is equal to the total ion- 
ization energies Va/b{°°) required to remove an electron 
from one atom and add it to the other. V then follows the 
Coulomb law until the atoms come quite close, when the 
atomic wavefunctions start to penetrate each other. At 
very small distances, the difference in energy is given by 
the energy level difference of the atomic wavefunctions, 
D = 1/2(Ea — Eb)- The transition from the Coulomb 
regime to the energy level difference regime may be ex- 
trapolated using a function of the following form: n 




±D-V A /b(°o 



- Va/b(c 



T-A/B 



ATre {±D-V A /B(™) 



(26a) 



(26b) 
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The length (1a/ b defines the length scale at which the 
transition between covalent and polar bonding occurs. 
Since this length is largely determined by the atomic size, 
it is typically close to the bond length. Still, the effect 
of Va/ b {oo) can often be ignored, since the effects of the 
energy level difference and the covalent bond combined 
tend to be stronger at atomic separations equal to the 
bond length. 

VII. CONCLUSIONS 

While many textbook explanations of chemical bonds 
try to visualize the bond only in its spacial distribu- 
tion, misleadingly suggesting that the source of the bond 
energy could be the slight increase in electron density 
between the atoms, the approach presented here clearly 
identifies the kinetic energy term as the dominant con- 
tribution to homopolar bond formation. In the spirit of 
the tight binding approximation, one can consider the 
atomic confinement of the electrons as much stronger 
than the interatomic interactions, and thereby arrive at 
a bonding potential by simply calculating the expecta- 
tion values of the kinetic energy for electron pairs evenly 
distributed between the atoms. This type of bonding 
potential, which uses only the unchanged atomic wave- 
function and the fundamental kinetic energy term of free 
electrons, not only gives the right order of magnitude for 
the bond energy, but also reproduces the spatial poten- 
tial of the bond, with realistic results for bond length and 
bond elasticity. 

Although a quantitatively accurate calculation may 
only be achieved by including a more detailed descrip- 
tion of polar effects than discussed here, it should be 
pointed out once more, that purely kinetic matrix ele- 
ments do give highly accurate results in semiconductors 
and other simple crystals. Furthermore, the free elec- 
tron behaviour of electrons in metals and the relation 
between covalent and metallic bonds may be understood 
better in the light of these considerations. In fact, the 
major difference between the two type of bonds is not the 
mechanism of bonding itself, but rather concerns the fact 
that metallic bonds are not directed as are bonds involv- 
ing p type wavefunctions. Instead they couple equally 
well to all neighbours, resulting in the strong derealiza- 
tion of electrons which makes metals conductors and in 
their relatively high plasticity compared with the brittle 
covalent bonds of e.g. semiconductors. 

In this manner, it is possible to find numerous con- 
nections between the very fundamental laws of quantum 
physics and the chemical and physical properties of the 
world surrounding us. Indeed, as technology advances it 
is important to remember that science not only tells us 
how to do things, but also why things are as they are, al- 
though this often teaches us more about limitations than 
about possibilities. Even the best artificial materials will 
not be orders of magnitudes removed from the typical 



properties dictated by constants such as Ti and e, just as 
no signal will ever travel faster than the speed of light. 
However, knowledge of our limitations may often proof 
more important than know-how, and in this sense, sim- 
plified models as the one presented here can proof to be 
extremely useful. 
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Table 1: 







bond 




bond 




bond 




Molecule 


Z 


length 


a 


energy 


a 


elasticity 


a 






in 1(T 10 to 


in 10- 10 m 


in kj/mol in 


10- 10 m 


in N/m 


in 10~ 10 m 


H 2 


2 


0.74 


0.23 


432 


0.49 


520 


0.36 




2 


2.67 


0.83 


110 


0.96 


130 


0.50 


Na 2 


2 


3.08 


0.96 


72 


1.19 


170 


0.47 


N 2 


6 


1.10 


0.34 


843 


0.60 


2260 


0.32 


2 


4 


1.20 


0.38 


494 


0.64 


1140 


0.35 


F 2 


2 


1.42 


0.44 


140 


0.85 


450 


0.37 


Cl 2 


2 


2.00 


0.63 


240 


0.65 


320 


0.40 


Br 2 


2 


2.28 


0.71 


190 


0.73 


240 


0.43 


CO 


4 


1.13 


0.35 


1071 


0.44 


187 


0.55 


NO 


4 


1.15 


0.36 


678 


0.55 


155 


0.57 


HCl 


2 


1.27 


0.40 


428 


0.49 


48 


0.65 



Table 1 shows the experimentally determined properties of some biatomic molecules, together with the value of the 
parameter which would reproduce this property using Gaussian wavefunctions. Note that despite the simplicity of 
the model, all values of a are of the same order of magnitude. 
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Figure 1: 

Figure la) shows the change in kinetic energy as a 
function of bond length for the bonding state. The inter- 
atomic distance d is given in units of 2a and the energy 
difference from the unbound state is given in units of 
h/4ma 2 . 

Figure lb) shows the corresponding k-space distribu- 
tions |^(fc) 2 | as a contour plot, the momentum Uk is 
given in units of Ti/a. 

The minimum in bond energy occurs when the k-space 
distribution is most narrow, just before the side maxima 
appear. 

Figure 2: 

Figure 2a) shows the change in kinetic energy as a 
function of bond length for the anti-bonding state, all 
units are as in Figure 1. 

For the anti-bonding state, the k-space distribution 
shows two peaks at small distances d. It is therefore 
much wider than the bonding state distribution. As d 
increases, it narrows, even though side maxima appear. 

Figure 3: 

Figure 3 shows the bonding potential of the ppa bond 
calculated for the p symmetry Gaussian wavefunctions in 
the xy plane. The unit of length is 2a. 
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Figure 1 : 



a) bonding potential 
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b) k-space distribution 
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a) anti-bonding potential 
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b) k-space distribution 
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b) three dimensional representation 



